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The ability to track very weak signals from distant spacecraft is limited by the phase 
instabilities of the received carrier signal and of the local oscillator employed by the 
receiver. These instabilities ultimately limit the minimum loop bandwidth that can be 
used in a phase-coherent receiver , and hence limit the ratio of received carrier power to 
noise spectral density which can be tracked phase-coherently. This article presents a 
method for near-real time estimation of the received carrier phase and additive noise spec- 
trum , and - optimization of the phase locked loop bandwidth. The method was used with 
the breadboard DSN Advanced Receiver to optimize tracking of very weak signals from 
the Pioneer 10 spacecraft , which is now more distant than the edge of the solar system. 
Tracking with bandwidths of 0.1 Hz to 1.0 Hz reduces tracking signal threshold and 
increases carrier loop SNR by 5 dB to 15 dB compared to the 3 Hz bandwidth of the 
receivers now used operationally in the DSN. This will enable the DSN to track Pioneer 10 
until its power source fails near the end of the century. 


I. Introduction 

The ability to track very weak signals from distant space- 
craft is limited by the phase instabilities of the received car- 
rier signal and of the local oscillator employed by the receiver. 
These instabilities ultimately limit the minimum loop band- 
width that can be employed in a phase-coherent receiver, 
and hence limit the ratio of received carrier power to noise 
spectral density which can be tracked phase-coherently. To 
achieve the best performance in a telemetry system, the spec- 
tral density of the phase noise should be known, and the pa- 
rameters of the phase locked loop should be optimized. 

At the present time, the Deep Space Network (DSN) uses 
two distinct operational telemetry receivers. These receivers, 
designated Block III and Block IV, respectively, employ sec- 


ond order analog carrier tracking loops, with a minimum loop 
bandwidth of 3 Hz. A 1 Hz bandwidth is possible with the 
Block IV receiver when tracking 2.3 GHz (S-band) signals, but 
this bandwidth is seldom used due to receiver oscillator noise, 
and because of the relatively large phase error caused by the 
inability of narrow bandwidth second order loops to track 
dynamics introduced by the rotation of the earth. 

One of the important advantages of digital phase locked 
loops over their analog counterparts is in the flexibility and 
accuracy of setting and modifying loop parameters. In the 
DSN Advanced Receiver, which is currently in breadboard 
form, the loop filter characteristics can be altered in real time. 
This capability motivates the search for an adaptive algorithm 
to optimize loop parameters, particularly loop bandwidth, in 
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real time. Optimization is especially important when tracking 
weak signals generated by noisy oscillators, such as the signal 
currently being received from Pioneer 10. When the loop is not 
stressed by dynamics, or equivalently when the loop order is 
high enough to track dynamics with negligible phase error, 
then the optimum loop bandwidth may be defined as that 
bandwidth which minimizes the steady state rms phase error 
within the loop. At its optimum value, the loop bandwidth is 
typically wide enough to track out phase jitter in the received 
carrier, yet narrow enough to effectively limit the amount of 
thermal noise allowed into the loop. By reducing the carrier 
phase jitter, radio loss is reduced and data recovery improved. 

In order to optimize loop bandwidth, it is first necessary to 
model the phase noise properly. The significant components of 
phase noise are due to the spacecraft transmitter, propagation 
effects, and VCO noise at the receiver. Spacecraft transmitter 
noise depends on the mode of transmission, that is, one-way, 
two-way, or three-way. With one-way transmission, a free run- 
ning oscillator on the spacecraft determines the transmitted 
phase. With two- and three-way transmission, the spacecraft 
phase locks to an uplink signal, multiplies its frequency and 
phase by a fixed ratio, and uses the resultant as the downlink 
carrier. If the downlink receiver and uplink transmitter are at 
the same station, then the mode is termed two-way, while the 
term “three-way transmission” is reserved for the case when 
the downlink receiver is at a different station. In the absence 
of propagation effects, one-way reception is usually dominated 
by spacecraft oscillator noise, whereas uplink transmitter noise 
normally dominates two- and three-way reception. A signifi- 
cant propagation effect is that due to solar scintillation, caused 
by charged particles in the vicinity of the Sun. This effect 
dominates when the spacecraft is more distant than the Sun, 
and the Sun-earth-probe (SEP) angle is small. Clearly, all of the 
above effects depend on a variety of factors; hence any real 
time characterization of the phase process can be of great 
value during reception. In this article we present a method for 
estimating the total phase spectral density in a phase locked 
receiver, and specifically in the DSN Advanced Receiver. The 
technique involves measuring the spectral density at the phase 
detector output, and extracting the individual spectral compo- 
nents by processing the estimated received spectrum. Spectral 
estimates are obtained by the application of fast Fourier trans- 
forms (FFTs) to the phase detector output. Given knowledge 
of the closed loop transfer function, the phase spectrum, 
receiver noise spectrum, and received carrier power can all be 
estimated. The optimum loop bandwidth is determined from 
these estimated components. Analytic expressions are derived 
for the case when the phase noise spectral density varies as 
/" a , a model often quoted in the literature. With a equal to 
three, the model can be used to describe oscillator phase 
spectra at low frequencies, while for solar scintillation spectra 
a value of a = 8/3 is appropriate. 


Finally, we present measurements of phase spectra, and 
determine optimum loop bandwidths for Pioneer 10 telem- 
etry. Data for both one-way and three-way transmissions are 
included. It is shown that the Advanced Receiver can reduce 
the minimum tracking threshold for Pioneer 10 by as much 
as 5 dB to 15 dB compared to the tracking threshold of 
current operational receivers. This is accomplished by using 
0.1 Hz to 1.0 Hz loop bandwidths, compared to the 3 Hz 
bandwidths currently used operationally. 


II. Mathematical Models 

In the following development, we assume that the received 
signal is of the form 

s(f) = \Jl A cos(a> 0 f + 0(f) + 0) 

where A is the signal amplitude, co 0 is the carrier radian fre- 
quency, 0(f) is a random phase process, and 0 is a random 
initial phase. Multiplication by \Jl has the advantage that the 
total carrier power becomes P c - A 2 , and the slope of the 
phase detector S-curve near the origin becomes A . This signal 
representation also allows direct comparison between the 
baseband phase spectral densities within the phase locked 
loop and their RF counterparts, where oscillator phase spectra 
are generally measured. Since in the current application the 
signal amplitude remains essentially constant for long periods 
of time (and hence can be estimated with great accuracy), we 
assume that A is known. 

A. Derivation of Observation and Phase 
Error Spectra 

A linear model of the carrier loop is shown in Fig. 1. This 
model is accurate when the loop is locked, and operating with 
small rms phase errors. Under these conditions, the loop is 
continuously generating estimates of the received phase 
process 0(f ), denoted 0(f), and subtracting the estimate from 
the received phase. The resulting phase error process 0(f) = 
0(f) - 0(f) is multiplied by the slope of the phase-detector 
S-curve at zero frequency (in our case that value is “ A ”), and 
filtered by the loop filter in the presence of additive receiver 
noise n{t). The output of the loop filter is an error process 
e(t) which is used to control the frequency of the voltage- 
controlled oscillator (VCO), or the equivalent numerically 
controlled oscillator (NCO) typically employed in digital loops. 
Instabilities within the receiver VCO are modeled as an addi- 
tive phase process 0 r (f) that is added to the VCO output. At 
this point it is useful to introduce the Heaviside differential 
operator p ~ djdt, which allows relating the output “z” of a 
linear filter to its input “x” in the time domain as z(f) = 
F(p) x (f), where F(p) is a ratio of polynomials in p . The 
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transfer function of the filter in the frequency domain 
becomes F(jaS) = F(p)\ p=fo ,. (Subsequently we shall use 
both radian frequency cc and frequency /= co/27r, as appro- 
priate.) Suppressing time dependence and making use of the 
Heaviside differential operator, the loop equations can be 
expressed as 


<t> = e -e 

(la) 

y = A<j> + n 

(lb) 

-HTI' 

(lc) 

0 = x + \p r 

(Id) 


where K is the filter gain, G the VCO gain, and y the noise- 
corrupted observable. In general, the received phase consists 
of phase processes due to modulation, doppler, and transmitter 
VCO instabilities. Any initial phase offsets are assumed to have 
been tracked out by the receiver. For the type of signals typi- 
cally observed by the Advanced Receiver, the use of sub car- 
riers can be assumed, which implies that modulation sidebands 
are far removed from the carrier and hence may be ignored. 
Therefore, the received phase can be modeled as 0 = d + \p t 
where d is due to possible doppler rates and 0 r describes 
the effects of additive transmitter VCO instabilities. Defining 
the transceiver phase process A\p = \p t - \j/ r and letting B(p) 
- KGF(p)/p , Eqs. (la), (lc), and (Id) can be used to obtain 

6 = B(p)y+\p r = AB(p)tp + B(p)n + \p r (2a) 


2 , W 

S 0 (f) = II - H(j2nf)\ 2 S ( (f) + I H(;2 nf) | 2 (5a) 

S (/) S (/) 

-Ar = II - HU2nf) I 2 SJf) + 1 1 -H{j2«f) | 2 
A 2 S A 2 

(5b) 

where f = d + A0. If the doppler and transceiver phase 
processes are independent, then S f (/) = S d (f) + 

Note that in Eq. (5b) we divided both sides of the equation 
by A 2 in order to emphasize the similarities between the 
observation and phase error spectra, when both are expressed 
in the same units. Aside from a scaling factor, the only differ- 
ence between the two equations is that the additive noise 
component in the phase error process is filtered by H(j2irf), 
whereas in the observation process it is effectively filtered by 
(1 -tf(/27T/)). 

Now let us consider the received phase process, specifically 
in the case where d is small. The received phase spectrum then 
depends on the transmitted phase spectrum and on the proper- 
ties of the channel. The free space channel is generally con- 
sidered to be an additive white Gaussian noise channel which 
does not introduce phase distortion. For now we will neglect 
the effects of the sun and the solar corona on the received 
phase spectrum, as these will be discussed in a later section. 
Therefore, in this case the phase stability of the received car- 
rier is dominated by the phase stability of the transmitted 
carrier. Two transmission modes commonly used in the DSN 
will be examined in this analysis. 


or 


0 = d + A 0 -AB(p) 0 - B(p) n (2b) 


0 = 


(1 +AB(p)\ 


(d+ Ai/0- 


( AB(p) 1 n 
(1 +AB(p)\ A 


(3a) 


In the one-way mode, the spacecraft oscillator is the source 
of the transmitted carrier. The term one-way stems from the 
notion of the frequency reference traveling only one way, 
from the spacecraft to the ground station. In the one-way 
mode, the phase noise spectrum is normally dominated by a 
term varying as inverse f 3 , for frequencies very close to the 
carrier. 


y = A 


1 +AB(p) 


(d + AiP) + 


1 

1 +AB(p) 


(3b) 


Since the closed-loop transfer function of the loop is 


HO w) 


AB(j go) 

(1 +AB(/u > )) 


(4) 


the power spectral densities of the observation and phase 
error processes can be represented as 


The other transmission mode discussed here is known as 
three-way transmission. In this mode, a carrier reference is 
transmitted to the spacecraft from a ground station. The 
spacecraft receiver phase locks to this signal, multiplies the 
frequency by a known ratio, and uses the resultant as the 
downlink frequency. Thus the received phase depends on the 
phase characteristics of the uplink oscillator and of the space- 
craft phase locked loop. In our experiments with the Pioneer 
10 spacecraft, transmission in the two-way mode is not feasi- 
ble since the round trip light time is roughly 1 1 hours 30 min- 
utes, so that by the time the transmitted reference is returned 
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by the spacecraft, the transmitting antenna is on the other 
side of the world. 

For three-way telemetry, the received process spectrum 
cannot be described by a simple inverse / 3 model, in general. 
Now the effects of the transmitter tracking loop on the uplink 
spectrum must be properly taken into account. A linear model 
of the tracking loop on the spacecraft appears in Fig. 2. This 
model is similar to the one describing the receiver tracking 
loop, except that now the loop’s output is its estimate of the 
uplink phase. The subscript “f” denotes “transmitter,” dis- 
tinguishing the components of this loop from those of the 
receiver. The transmitter loop tracks the uplink phase, using it 
as the reference for the downlink telemetry. Thus, the phase 
of the downlink carrier is identical to the loop’s estimate of 
the uplink phase. The power spectral density of the down- 
link transmitted phase process B t can be derived with the same 
techniques used to derive the spectral density of the receiver 
phase process. The result is 


SAf) = \H(j2nf)\ 2 

t 



+ \l-H t U2nf)\ 2 S^if) (5c) 

where 6 U is the random phase of the uplink oscillator multi- 
plied by the spacecraft transponder ratio. Note that the uplink 
phase process and the internal noise are both filtered by the 
closed loop transfer function, while the transmitter’s phase 
process is effectively filtered by the complementary filter 
function (1 - //(/27r/)). In our case Sq u — 10“ but 

1 1 - H t ( j2nf) | 2 < 10 -4 at the low frequencies of interest, so 
that typically Sq u is the dominant component of the downlink 
phase process. 

The transmission channel may also affect the spectral 
purity of the carrier phase. A good example of a prominent 
channel effect is provided by the solar corona [1] . The result- 
ing phase degradation can be appreciable, often negating the 
benefits of highly stable earth-based oscillators operating in 
three-way mode. The power spectrum of the phase degrada- 
tion for this solar effect is typically inversely proportional to 
the 8/3 power of frequency. 

B. Closed Loop Transfer Function and Its Effects 

Next we examine the structure of the closed-loop transfer 
function in greater detail, and consider its effects on various 
phase process and additive noise models. 

The Advanced Receiver employs second and third order 
loops. Approximating the rapidly sampled data loop by a 


continuous time loop, the transfer functions of the loop 
filter and of the resulting closed loop for a third order, type 
III loop are [2] 


1 + t s 

F(s) = L + 


V T 2 T 3 r 


(6a) 


r(k + t s+(t s) 2 ) 

H(s ) = (6b) 

r (£ + t 2 s + (t 2 s) 2 )+(t 2 s ) 3 

where r, , r 2 ,r 3 are constants, r = AIct%It 1 , k = r 2 /r 3 and s 
is the Laplace transform variable. The steady -state transfer 
function is obtained by letting s - j co in Eq. (6b), yielding 


H{i w) 


1 - H{j co) 



(7a) 


(7b) 


where 


* = B l h 


r (r-k+ 1\ 
4 V r-k ! 


( 8 ) 


These transfer functions can also be expressed in terms of the 
normalized radian frequency co = co/Z? L , indicating that B L 
simply scales the frequency variable. Note that the loop re- 
duces to a second order type II loop when k = 0. The squared 
magnitudes of the transfer functions defined in Eq. (7) are of 
fundamental importance in determining loop performance. 
These functions are 


^( S ) - \h(ju)\1 =b v - r2(k < ° 2 - )+S52( — 


' L w r 2 (k - £j 2 g 2 ) 2 + (o5rg - oj 3 g 3 ) 2 

(9a) 


FJu 3) - |1 -H(j u)\ 


2 _ = CoV 

r 2 (ic-Z5 2 g 2 ) 2 +(S$rg-(Z 3 g 3 ) 2 
(9b) 
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Plots of these functions are shown in Fig. 3 for r - 4, and 
k - 0.0 and 0.25. Note that for k < 0.5 and r > 2,F 2 (o5) 
is a monotone increasing function of <3, which implies that 
B l 2 > B lx , F 2 (co/i? L1 ) > F 2 (o;/Z? L 2 ). This property can be 
used to establish the existence of an optimum loop bandwidth. 

C. Oscillator Phase Noise 


trum always occurs at a frequency that is a linear function of 
the loop bandwidth. This can be established by differentiating 
| 1 - H{ju) | 2 Sf(f) with respect to c o and setting the result 
equal to zero. The maximum occurs at frequency / = /*, where 

f* = C 0 (r,k)B L (13) 


Oscillator phase noise can be modeled as [3] 


S*(f) 





/ 4 


\f\<F u 


( 10 ) 


\f\>F u 


where F u is an upper frequency cutoff and the S t , i = 1,4, 
are constants. For our interest, we observe that if the trans- 
mitter and receiver oscillators are running independently, 
then their spectral densities add, yielding • = S$ f (f) 

+ Each component is typically of the form given in 

Eq. (10), although the corresponding constants may be dif- 
ferent for the two oscillators. The low-frequency behavior 
of an oscillator is usually dominated by “frequency flicker 
noise” [4] , giving rise to the well known inverse / 3 behavior 
of phase fluctuations. Assuming that in the region of interest 
frequency flicker noise dominates, and further assuming 
negligible doppler contribution, the power spectral density of 
the phase process for one-way telemetry can be modeled as 


(the derivation is shown in the Appendix). For (r = 4, k - 0.25), 
/* = 0.127 B l ; for a second order loop with parameters (r = 4, 
k ~ 0.0), /* = 0.147 B l . This result can be used to establish 
whether or not an observed phase process is adequately 
described by the/’ 3 model. 

D. Optimization of Tracking Loop Bandwidth 

The performance of phase-locked loops can be assessed in 
terms of the total rms phase error a ^ present in the loop. The 
total mean-squared phase error is the integral of the phase 
error spectral density defined in Eq. (5a): 

= 2 SJf) df ± aj + a* (14) 

The components of the total phase error oj and a 2 are 
defined as 

o) = 2 J |1 -HU 2 irf) | 2 S t (f)df (15a) 
and for white noise spectra of the formS n (f) = N 0 /2, 


S, 

Sy(f) - — , (11) 

S I/I 3 

This model may also be applied to three-way telemetry, if 
the spacecraft loop bandwidth is so great that the uplink phase 
process term dominates. 

At low frequencies the additive receiver noise is generally 
white. Denoting its single-sided spectral level by N 0 yields the 
representation 

N 

s n (f ) as ~Y (12) 

With these approximations, the process, additive noise, and ob- 
servation spectral densities appear as in Fig. 4, for a third order 
loop with parameters as before. 

Note that for filter functions of the form defined in Eq. (7) 
and inverse / 3 phase processes, the peak of the process spec- 


(15b) 

Note that for white noise the spectral level is constant; hence 
the component of the variance due to white noise increases 
linearly with loop bandwidth. The process component depends 
on the power spectrum of the phase, which is generally not 
constant. However, since the power spectral density is a non- 
negative function of frequency, it follows that for B L2 >B L1 , 
S t ( co) F 2 (cojB L1 ) > Sf( go) F 2 (oj/B L2 ). Since this inequality 
also holds for the integral, the component of the variance due 
to process noise is a decreasing function of loop bandwidth, 
for all valid power spectral density functions. It follows from 
the above argument that a loop bandwidth exists which mini- 
mizes the total mean squared phase error. 

Loop bandwidth optimization can be achieved for arbitrary 
phase spectra by computing estimates of the total phase error 
as a function of loop bandwidth, and selecting the bandwidth 
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corresponding to the minimum of this function. For the case 
of additive white receiver noise, this requires knowledge of the 
parameter N 0 /2A 2 f and of the function |1 -H(j 2rrf)\ 2 £$.(/), 
which is the integrand in Eq. (15a). If the power spectral den- 
sity of the observable and the density of the signal level A are 
known, then the required integrand can be determined from 
Eq- (5b): 


| 1 - H(j 2n f) | 2 SJf) = - 1 1 -H(j 2n f) l 2 (-^-J 

r A 2 \2A 2 J 

(16) 

A major simplification in the optimization algorithm results if 
the phase process spectrum is assumed to be of the form /~ a . 
This assumption is normally valid for one-way transmission, 
and may be accurate for three-way transmission as well, pro- 
vided that the instability is dominated either by uplink oscilla- 
tor or solar scintillation phase noise. Introducing the change 
of variables / = / /B L , it follows that 

<7 2 = 2 —> | 1 - H(j 2tt /) | 2 df 

= 2S 1 T 0 (a)£i 1 - Q) (17) 

where 7 0 (a) is the value of the integral. Using numerical inte- 
gration, we found that 7 0 (3) = 9.08 and 7 0 (8/3) = 5.87 for a 
third order loop with parameters r = 4, k = 0.25 and a loop 
bandwidth of 1 Hz. 

For the model, the minimum of the total phase error 
variance can be found by differentiating a £ with respect to B L , 
setting the result equal to zero, and solving for the optimum 
loop bandwidth 

3a 2 /N.\ 

d -^=-2(a-l =0 (18a) 


for loop bandwidth optimization requires only knowledge of 
oj, a 2 , and the bandwidth at which the measurements were 
made. Measurement of aj also allows determination of the 
phase process constant S l , via the equation 


*5 


l 


2 7 0 (a) 


(19) 


In a practical system, the required parameters and spectral 
densities are not known and hence must be estimated. The 
accuracy of the predictions ultimately depends on the accu- 
racy of the estimates. The choice of best estimator structure 
often depends on the application. The selection of the esti- 
mator structure was guided by the fact that in our application 
the amount of available samples far exceeded the requirements 
imposed by resolution constraints. 


E. Power Spectrum Estimator Using Bartlett’s 
Procedure 

The power spectral density estimates in this article use 
Bartlett’s technique of averaged periodograms, which involves 
a trade-off between smoothing and spectral resolution. Bart- 
lett’s approach can be used to advantage whenever the avail- 
able record length is so great that the required spectral reso- 
lution can be achieved with a small fraction of the available 
samples. Thus, if the total record length is K but N « K 
samples satisfy the resolution requirements, then M = K/N 
periodograms may be averaged to obtain a smoothed estimate 
of the spectral density. Our interest in periodograms stems 
from the fact that FFTs can be used to compute sampled 
versions of the desired periodograms quickly and efficiently. 


The periodogram I N ( co) associated with a sequence y{ri) 
of length TV is defined as 

I Jo,) = L\Y(e^)\ 2 (20a) 


where 


yields 

1 /<* 

( (a - l)a 2 ) 1/01 


(18b) 

Figure 5 shows the behavior of a? and o 2 as a function of loop 
bandwidth, as well as their sum, for a = 3 and parameter values 
S x = 477 X 10~ 4 and (7V 0 /2^4 2 ) = 2 X 10“ 2 . This algorithm 



N - 1 

Y(e i(JJ ) = £ y(n)e~ iu}n w = 2 irf (20b) 

n~0 


is the discrete Fourier transform of the sequence y(n). The use 
of periodograms in spectral estimation can be justified on the 
grounds that in the limit as N approaches infinity, the expect- 
ed value of the periodogram approaches the power spectral 
density of i.e., 

lim ^[/^(w)] = S (oj) (21) 

N-»°c 
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Thus, in the limit the expected value of the periodogram is the 
desired power spectral density. 

The variance of the periodogram as a function of frequency 
provides an indication of the amount of random variation 
about the mean of the estimate. For a real Gaussian sequence 
with power spectral density S y (<jo), the covariance is approxi- 
mately [5] 



Note that the standard deviation of a single periodogram is at 
least as great as the spectral level itself, at any frequency. At 
frequencies away from zero, the standard deviation of I N (oo) is 
well approximated by I N ( co). Therefore, a single periodogram 
does not provide a very useful estimate of the power spectral 
density; hence averaging is usually required to reduce the vari- 
ance of the spectral estimate. 

The FFT can be employed to generate samples of I N ( go). 
With 


» £sj<*) H ♦V.J « vtj 

( 24 ) 

where bk ly k 2 * s Kroenecker delta (this function is non- 
zero only for k x = , where its value is one). In particular, the 

variance is 

var [£,(*)] = ~S 2 y (k) k* 0 (25) 

The spectral estimate samples S y (k) are independent random 
variables, by virtue of the fact that the final estimate is an 
average of M random variables. Thus for any k , averaging M 
independent periodograms reduces the variance of each esti- 
mate of spectral density by a factor of M relative to that of a 
single periodogram. 

F. Error Analysis 

The estimator described above cannot provide perfect esti- 
mates of the desired parameters by processing finite sequences 
(in fact, no estimator can). We determine the variance of the 

A 

estimates and B[, under the assumption that enough period- 
ograms were averaged to ensure that the errors in the under- 
lying parameters are but a small fraction of their actual values. 
Let dj and o'* denote the estimates of oj and o *, respectively. 
For small errors, we can write 

a\ = + e l = (26a) 


Y(k) = Y(e^)\ 


LO~2nkfN 


- °: + < = a:G+*J 


it follows that 


t N (k) / Ar( a; )l w = 27rJt/iV 


Since the samples y(ri) were assumed to form a Gaussian 
sequence, it follows that the FFT outputs y(/:) are Gaussian 
random variables. It is readily demonstrated that the random 
variables y(fc) are uncorrelated and hence independent by the 
Gaussian assumption. It follows that the periodogram samples 
I N (k) are independent as well. 

Let S y (k) denote the final spectral estimate, after averaging 
over M periodograms. It is evident that 


Assume that x ^ and x n are zero mean random variables, with 
rms values much less than one. Using Eqs. (18) and (19), the 
estimates of B£ and can be expressed in terms of their 
actual values as 


[l+x t ) 1/a 

B* = Bf{ T — L \ 

L L I 1 + x n 1 


lim /?[£ (*)] = lim i £>[/„ ,.(*)] = S Jk) 


s t = s,(l+* f ) 
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The second order statistics of S l can be obtained directly: |||. Pioneer 10 Results 


E[S 1 ]=S l (28a) 

var(a* ) 

var[Sj = S ] var(x ? ) = S] ■ — — (28b) 

(In fact, these second order statistics for S x are valid in gen- 
eral.) To obtain the statistics of expand the ratio contain- 
ing the error terms and ignore all but the linear terms: 

1 + *..) 1/a ( x 2 ) 1/a 

— j = | ( i + *$-)o + ~ - " j 

- ji +(* f -* M )J l/a 

( 29) 


Thus, for suitably small errors, we obtain the approximation 


This section presents the results of experiments in tracking 
the Pioneer 10 spacecraft with the breadboard DSN Advanced 
Receiver [6] . The data were obtained using the 64 meter an- 
tenna at DSS-14 in Goldstone, California. The data include 
both one-way and three-way transmission modes obtained at 
various SEP angles, illustrating the effect of the solar corona 
on carrier phase. Spectra were measured and loop bandwidths 
were optimized for the various conditions. 

The Advanced Receiver is a hybrid analog/digital receiver 
which makes digital estimates of the carrier phase and fillers 
these estimates with a digital second or third order tracking 
filter. The filter output drives an analog frequency synthesizer, 
closing the loop. The phase detector output provides the data- 
stream which is subsequently analyzed to determine the re- 
ceived spectral densities. The power spectral density at the 
phase detector output, S y , was obtained as described in Sec- 
tion II.E. These estimates were obtained in near-real time. 
Phase measurements as a function of time were clocked out 
of the Advanced Receiver to a control computer, where these 
measurements were recorded on floppy discs. After sufficient 
data were recorded, the disc was processed on an IBM PC. The 
processing was typically completed in less than five minutes, 
during which time data for the next tracking case was col- 
lected by the Advanced Receiver. 


s t - B t | < 30) 

It follows that 

E[B*\ = B* (31a) 

var t£*l = {^1 vaifcf-*,) 

(B* l \ 2 ( var(a^) var(o^) 

\/var(a!) v/var (?) I 

+ 2 3 } 

° 2 ) 

r " (31b) 


The upper bound in Eq. (31b) is useful when the correlation 
between x^ and x n is not well known. Actual numerical values 
will be determined in the following section, where the tech- 
niques for obtaining the parameters from the spectral esti- 
mates are examined in greater detail. 


A. Spectral Estimation Results 

Estimates of the power spectral density ofj>(w), S (k), were 
obtained by means of averaged periodograms, using the FFT. 
If the sampling rate satisfies the Nyquist criterion, the result is 
also a good estimate of the spectrum of the continuous process 
y(t). Results for one-way transmission are shown in Figs. 6a, b, 
and c. The phase spectral density estimates are represented by 
the solid curve (the points of the FFT output are connected 
by straight lines). For the first case, Fig. 6a, the number of 
points in the sequence was N = 128 and M = 25 records were 
averaged. The sampling rate was 8 samples per second. Since 
for real samples the positive and negative portions of the 
spectrum are identical, only the positive portion is displayed. 
For purposes of interpolation, 128 zeros were appended to 
each record prior to processing, and a 256 point FFT per- 
formed (this does not change the total number of independent 
samples in the FFT output, which remains N). In Fig. 6, a 
third order loop with parameters r = 4, k = 0.25 was used, with 
bandwidths of 0.5 Hz, 0.8 Hz, and 2.0 Hz (the 0.8 Hz band- 
width is nearly optimum). The contribution of the process 
spectrum near the origin dominates for the 0.5 Hz loop of Fig. 
6a; it is still present in the 0.8 Hz loop (Fig. 6b), but is not 
discernible for the 2.0 Hz loop (Fig. 6c). The dominance of 
the white noise component is apparent in all cases at frequen- 
cies greater than the loop bandwidth. 
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The above observation suggests the following procedure for 
estimating the desired components of the spectral density: 

(1) Estimate the normalized spectral level NJ2A 2 by aver- 
aging the discrete points of the spectral estimate over 
the upper half of the spectrum. (Henceforth, cap and 
Est [* ] shall both denote “estimate.”) 

(2) The function |1 -//( jk) |' 2 S n (k) is estimated, assuming 
the filter transfer function is known. 


shown in Fig. 7, again around the optimum loop bandwidth, 
which in this case was roughly 0.6 Hz. For this day, the SEP 
angle was 11.8 degrees. This resulted in process noise due to 
solar scintillation which was almost as high as the one-way 
phase noise due to the spacecraft oscillator. Estimates of the 
various system and model parameters can be obtained from 
these fundamental estimates. The following examples serve to 
illustrate some techniques for deriving the desired parameter 
estimates. 


(3) An estimate of the filtered process spectrum is ob- 
tained from Eq. (5b) as 


Est j |1 - H{jk) | 2 S f {k) j = 1 j S y (k) - |1 - H(jk) 1 2 S n (k ) j 

(32) 

(4) Finally, having decomposed the estimated observation 
spectral density into its components, estimates of the 
underlying spectral densities are obtained. 

Because the signal amplitude A is assumed known, it fol- 
lows that 


S n (k) = A 2 



(33) 


Since a 27V-point FFT was performed on an Appoint sequence, 
the resulting A^/2-point average has N/4 degrees of freedom 
and hence may be treated as the average of N/4 independent 
random variables. Taking into account the M records that were 
averaged to obtain the spectral density estimate, the mean and 
variance of this estimate are 



The smooth dashed curves in Figs. 6a through 6c are the esti- 
mates of receiver noise spectra filtered by (1 - H(jk)), while 
the rough dashed curves are the estimates of the filtered phase 
process spectral density. 

The same techniques can be applied to the case of three- 
way telemetry. Results for three-way tracking on DOY 141 are 


B. Phase Error Variance Estimates 

The components of the phase error variance can be found 
from the measured spectra. If the loop bandwidth is known, 
then from Eq. (15b) the variance of the additive noise compo- 
nent can be found since it is proportional to the normalized 
noise spectral level estimate: 


o 


2 

n 


= 2 B t 



(36a) 


Using the estimated spectral level, the variance of o 2 is approx- 
imately 


s m (2B F 


\ia 2 J ~ NM ^ n 


) 2 (36b) 


Computation of the variance due to the phase process 
requires integrating the process spectral estimate. For suitably 
great N , integration over the lower half of the positive spec- 
trum can be approximated by the sum of independent random 
variables as 


= 2 



II -//(/ 2tt /) | 2 (/)<// 


N 

-~2 
2 


= 2{A/ £ \1-H{jk)\ 2 fyk)} (37a) 

*= 0 , 2 , 4 , 


where F is the total frequency range and Af = F/N is the effec- 
tive spectral resolution of the FFT. Using Eqs. (32) and (35) in 
Eq. (37a), and since the spectral level estimates obtained from 
the upper half of the spectrum are independent of the lower 
half,, the variance of the estimate can be expressed as 
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var(a 2 ) s 4(A/) 2 




1 -H(jk)\‘ 


(37b) 


where we used the fact that the 7V/2*point estimate has only 
Nj 4 degrees of freedom. Near the optimum loop bandwidth 
the observed spectral level is essentially constant and equal to 
the normalized spectral level of the additive noise (see Figs. 6 
and 7). Upper bounding 1 1 -H(jk ) | 2 by one, Eq. (37b) can be 
approximated as 


a question mark. These estimates disagree with the other esti- 
mates by significantly more than the estimation error. These 
estimates depend on oj, which is the difference of two rela- 
tively large numbers whose difference is small. Thus the 
actual errors may be dominated by small biases that were not 
accounted for in the estimation of the larger quantities, one 
of which is (N 0 /2A 2 ). 

In Table 1 , the estimate 


P 

c 




(39) 


is displayed (in dB-Hz) instead of the normalized noise spectral 
level, in order to facilitate comparison with standard experi- 
mental measurements. 


var(o 2 ) 


=* n 2 

M 




= 4(A/) 2 

M 




(37c) 


For the Pioneer 10 data obtained on DOY 119 and 141 , the 
total phase error, i.e., the sum of the two components, is 
plotted as functions of loop bandwidth in Fig. 8, along with 
error bars of one standard deviation. Both one-way and three- 
way results are shown. It is apparent that in this case three- 
way reception is dominated by solar scintillation effects. 


The variances of parameter estimates are displayed in 
Table 2. These variances were computed using the constant 
spectral level assumption, and using the parameter estimates 
instead of their jtrue (but unknown) values in Eqs. (28b) and 
(31b). Because a 2 and are not highly correlated, we use the 
approximation var(a^) ~ var(a 2 ) + var(a^). The variances of 
Bl were computed by means of the upper bound in Eq. (31b) 
and hence tend to be too large. The standard deviation error 
bars in Figs. 8 and 9 were obtained directly from Table 2. For 
bandwidths smaller than the optimum, the actual standard 
deviations may be larger than indicated, because the spectral 
peak tends to be much greater than the spectral level of the 
additive noise used in the calculations. 


C. Optimum Loop Bandwidth 

Next we consider the estimation of the optimum loop 
bandwidth by means of Eq. (18): 

"*|l joe 

Bl = ^ B l (38) 

The estimate of the optimum loop bandwidth depends on the 
ratio of phase process and white noise variances, as well as on 
the loop bandwidth at which these variances are measured. 
Table 1 presents the estimated parameters for all data taken 
during Pioneer 10 tracking. At least three different loop band- 
widths were used each day. 

The estimates of optimum loop bandwidth are consistent 
on each day, except for the cases in Table 1 that are flagged by 


D. Phase Process Spectral Estimates 

Estimates of the phase process spectral density at unit 
frequency S x are presented for all of the Pioneer data in 
Fig. 9 and in Table 1. The spectral density at 1 Hz,Sj, was 
estimated by means of Eq. (19). The /~ 3 model was used for 
all one-way data, and the /~ 8 / 3 model was used for all three- 
way data. The SEP angle is indicated in Fig. 9 for each day 
data was collected. The results show that one-way data have 
spectral levels of S x = 1 to 2 X 10 -3 r 2 /Hz. 

For the three-way data, the results depend on SEP angle, 
indicating that here solar scintillation effects dominate. The 
measured spectral levels were approximately 10' 4 r 2 /Hz on 
DOY 128, when the SEP angle was 23.7 deg, and 10” 3 r 2 / Hz 
on DOY 141 at an SEP angle of 11.8 deg. Armstrong, Woo, 
and Estabrook [1] have measured solar scintillation dominated 
carrier spectra using the Viking spacecraft. At / = 0.001 Hz, 
they observed a spectral level of 10 3 to 10 s r 2 /Hz at 23.7 deg 
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and 10 3 4 to 10 6 * r 2 j Hz at 1 1.8 deg for one-way paths at a car- 
rier frequency of 2.3 GHz. Multiplying by 10 8 to convert to 
S j using the a - 8/3 assumption, and also multiplying by 2 to 
convert to three-way reception, the corresponding range for 
iSj becomes 2 X 10“ 5 to 2 X ICC 3 at 23.7 deg and 2 X 10" 4 to 
2 X 10“ 2 at 11.8 deg. Although extrapolation to frequencies 
as high as 1 Hz may not be accurate, we observe that our 
results are consistent with the previously measured values. The 
variation in the S j estimates evident in the three-way mode 
(Fig. 9) could be the result of actual solar scintillation, or 
could be due to modeling errors. For example, during three- 
way reception, the phase noise is due to a combination of 
spacecraft transmitter phase effects and solar scintillation; 
hence the effective exponent a may not be exactly 8/3, as was 
assumed in the calculations. In addition, the loop bandwidth 
in the Advanced Receiver is not precisely known. The esti- 
mates of S j are sensitive to both of these effects. However, 
since loop bandwidth optimization is less critically dependent 
on B l and a, our model is sufficiently accurate for the purpose 
of optimizing receiver loop bandwidth. Since the spectra 
measured in [1] tend to exhibit significant fluctuations, our 
method of measuring spectra in real time should prove to be 
very useful for optimizing loop bandwidths and for obtaining 
other useful system parameter estimates during telemetry. 


IV. Conclusions 

A method for optimizing loop bandwidths via spectral esti- 
mation for phase locked loop receivers has been developed and 
used to improve carrier tracking for the Pioneer 10 spacecraft. 
Estimates were made of the relevant spectral parameters and 
of the total noise power for one-way and three-way transmis- 
sion, and the variance of the phase error in the loop was mini- 
mized. Results obtained in the field were found to be in good 
agreement with theoretical values. In addition, by analyzing 
phase spectra, solar channel effects were detected and compen- 
sated for. Thus the Fourier analysis of the phase detector out- 
put was shown to be an important tool which in addition to 
bandwidth optimization can be used to monitor channel 
effects, transmitted carrier stability, and make estimates of 
relevant system parameters. 

It was specifically shown for Pioneer 10 that receiver loop 
bandwidths of 0.1 to 1.0 Hz are feasible, depending on the 
type of transmission (one-way or three-way), and depending 
on the SEP angle. When compared to the 3 Hz loop band- 
widths used with operational DSN receivers, it is clear that 
the narrower loop bandwidths achieve improvements of 5 to 
15 dB in carrier tracking threshold and carrier loop SNR. 
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Table 2. Variances of parameter estimates 
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Fig. 6. Spectral density estimates (one-way mode) 



Fig. 5. Phase error variance and components as a function 
of loop bandwidth (a = 3) 
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Fig. 8. Measured total phase error variances as a function 
of loop bandwidth 
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Fig. 9. Estimates of phase spectrum parameter 
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Appendix 

Maximum of Filtered Spectrum 


Here we demonstrate that for /~ 3 type phase processes the 
maximum of the filtered process spectrum |1 -//(/ 2n f)\ 2 
S^(f) is always at a given fraction of the loop bandwidth, for 
transfer functions of the form defined in Eq. (7) and inverse 
/ 3 process spectrum defined in Eq. (11). These functions can 
be expressed in terms of the radian frequency co = 27r/as 


= (2*) 3 — , (Al) 

I col 


and 


x 3 + ax + b = 0 



(A5) 


(A6a) 


(A6b) 


I 1 - //(/to) I 2 S f (w) = (277) 3 G(co-r,k,T 2 ) (A2) 


where 


Trigonometric solutions for x are obtained by letting x = 
mcos 6, with m - 2 yj-aj 3 . The solutions for are of the 
form 


G - cu 3 [co 6 + co 4 + a 2 co 2 + a 3 ] 


-l 


1 ^2 


r 2 (l-2fc) . . 

a, = — - — - (A3) 


r 2 k 2 

* 3 = T 


Setting 7 = co 2 , differentiating Eq. (A3) with respect to.y, and 
setting the result equal to zero yields the cubic equation 


y. = {Ivf C]{r,k)Bl (A7) 

Taking only the positive solutions, we obtain 

co* = 2tt C 0 (r,k)B L (A8) 

Direct evaluation yields the location of the maxima for r - 4 

k = 0: /* = C q (4,0)B l = 0A41 B l (A9a) 

k = 0.25: /* = C q (4,0.25)^ l = 0.127 ^ (A9b) 



Further, letting x =y + (c^/9) results in the simplified form 


These results are confirmed by the graphs shown in Fig. Al, 
where the location of the maxima are seen to be at the pre- 
dicted frequencies. This property helps to determine if the 
observed phase process is indeed an/ -3 type process. 
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Fig. A1. Filtered process spectra for second and third 
order loops 



